Spin Accumulation in the Electron Transport with Rashba Interaction 
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The non-equilibrium transportation of two-dimensional electrons through a narrow channel is 
investigated under the influence of the Rashba interaction. By introducing suitable lifetime in the 
Green's function, the average spin values can be calculated from the ballistic regime to the diffusive 
regime. It is shown that the spin accumulation is a combined effect of the spin current and disorders. 
In the diffusive regime, disorders offer a mechanism to stop the spin current and generate the spin 
accumulation. In the ballistic regime, spins are more spread out and do not have definite signs. 
Further consideration indicates that the inclusion of ferromagnetic spin-spin interaction increases 
the spin accumulation near the edge. 

PACS numbers: 72.25.-b, 72.25.Dc, 73.40.-c 



Generation and transfer of spins, preferably without 
using external magnetic field or magnetic materials, is 
an important issue in spintronics[l]. In semiconduc- 
tors, through the intrinsic spin-orbit interactions related 
to the Dresselhaus 2] or the Rashba^ mechanism, it 
is possible to generate, or even manipulate spins us- 
ing only an electric field. Recently, it is further pro- 
posed that the spin Hall current can be generated in bulk 
semiconductors |^ and semiconductor heterojunctionsjlj 
because of the spin-orbit interactions. Unlike an earlier 
proposed spin Hall effect involving impurities 6] , these 
two studies demonstrated that the spin Hall effect is pos- 
sible in pure samples without any impurities. Experimen- 
tally, the spin Hall effect has been observed in a number 
of systems pj . The key observation is the accumulation 
of opposite spin polarizations on two sides of the sample. 
While it is natural to associate the spin-accumulation 
with the spin Hall effect, for the intrinsic case, a satisfy- 
ing theory of spin accumulation still does not exist. 

In this paper, by introducing suitable lifetime in the 
Green's function, the spin accumulation will be calcu- 
lated from the ballistic regime to the diffusive regime. 
The origin of spin accumulation will be clarified. 

We start by considering a two-dimensional electron 
system confined in a channel 

H = -^x/ 2 +^z-[axp}+V c + V I . (1) 

Here a characterizes the Rashba iteraction, a are Pauli 
matrices and Vq is the hard-wall potential that confines 
the electron to — W/2 < y < W/2 . Vj is the potential 
due to point disorders, Vr(r) = J2i v ob~( r — Rz)0- We 
shall use 7TT to denote averaging over disorders and 
to denote the concentration of impurities. Since H pos- 
sesses the time-reversal symmetry, the Kramer's theorem 
implies that each eigenstate is doubly-degenerate. Physi- 
cally, the degeneracy represents two opposite propagating 



directions along x axis. Therefore, we shall denote the 
Kramer degeneracy by A = ±, representing propagating 
along ±x respectively. 

To find the spin-accumulation (S z (r)) (= s z (r)), we 
first note that because cr y H(x,y)a y = H(x,—y), for 
each eigenstate i/j n x(x,y), one obtains a v t/j n \(x,y) = 
ipn\(%> —y) and thus it follows that s z (r) is antisymmet- 
ric in y. This simple argument, however, can not as- 
sign a definite sign for s z (r). Numerics is then necessary 
to pin down the sign. It is important to note that for 
each energy, s z (r) have opposite signs for opposite prop- 
agating directions. As a result, spin-accumulation can 
only happen when the system is not in equilibrium, i.e., 
there must be a biased potential between two electrodes. 
Therefore, one has s z (r) =£„ « + (r) | Sf | ^ + (r) ) , 
where ip^ + (r) denotes eigenstates with A = + and the 
energy E n is restricted to the regime eV\ < E n < eV 2 
with V\ and V% being the electric potentials of electrodes. 
By inserting J^y 2 dES(E — E n ) and using the identity 
lm[l/(E - E n +ie)] = -iirS(E - E n ), we find 



(S,(t))= dE^—hx>Tr[S s G+(E,r,T)]j, (2) 

where G+ (E, r, r) is the retarded Green's function in 
the A = + channel and Tr is the trace over the 
spin space. In the Born approximation, G°L (E, r, r) = 

EJC+(r)]*<C(r)/[£ - E n + i 7 N + (E,r) + ie] with 
7 = tIjUq, <fi n+ (r) being the energy eigenstate and 
N + (E, r) being the local density of states in the absence 
of disorders^. Because k x is a good quantum number in 
the absence of disorders, A = + corresponds to positive 
k x . For a given k x , </>„+(r) is a linear combination of 
four waves with k y = where E = h 2 k±/2m ± ak± 

with k\ = ^ + (fc^) 2 . The hard- wall boundary condition 
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yields the following equation and further selects ky 
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Once are found, <j) n + and thus G°f (E,r,r) are ob- 
tained. Clearly, ky determine how s z (r) oscillates for 
each state and in general, there is no obvious spin ac- 
cumulation for each state. This is similar to the charge 
Hall effect in which there is no charge accumulation for 
single particle states and one needs Coulomb interaction 
to stop charges and generate charge accumulation. Fig. 
1 shows the numerical results of s 2 (r) for different 7. 
By setting 7^0, one reaches the ballistic limit where 
spins are more spread out, and furthermore, as shown 
in the inset of Fig. 1, depending on parameters, s z (r) 
may even switch the sign. Nonetheless, as one turns 
on disorders, s z (r) begins to accumulate near the edge 
and always switches to be negative on the left hand side 
(LHS) and positive on the right hand side (RHS). Obvi- 
ously, disorder is mechanism for generating spin accumu- 
lation. The reason can be found by exploring the con- 
tinuity equation ds z (r,t)/dt = —s z (r,t)/r — V • 3 z (r,t) 
where r is the effective diffusion time due to disorders 
and 3 Z is the spin current. In the steady state, one ob- 
tains s z (r) = — tV • J z (r), which implies that if J z flows 
from right to left, s z (r) is negative on LHS and positive 
on RHS, in consistent with numerical results. Thus the 
spin diffusion generates the spin accumulation that stops 
the spin current. 

In addition to the spin diffusion, the spin-spin inter- 
action, JJ2 cr i ' °j ' further enhances the spin accumula- 
tion. Fig. 2 shows our numerical results based on the 
self-consistent mean-field theory . Clearly, ferromagnetic 
coupling amplifies the oscillation of the spin density and 
increases the accumulation. For antiferromagnctic cou- 
pling, the oscillation has a shorter length scale because 
neighboring spins tend to be antiparallel. 

In conclusion, we have investigated the spin accumula- 
tion from the ballistic limit to the diffusive regime. The 
origin of the spin accumulation is clarified. This research 
was supported by NSC of Taiwan. 



[4] S. Murakami, N. Nagaosa, and S. C. Zhang, Science 01 
(2003), p. 1348 ; Phys. Rev. B 69 (2004), 235206. 

[5] Jairo Sinova, Dimitrie Culcer, Q. Niu, N. A. Sinitsyn, 
T. Jungwirth, and A. H. MacDonald, Phys. Rev. Lett. 
92(2004), 126603. 




FIG. 1: Spin accumulation for different 7. Here eVi = 50, 
eVi = 0, and koW = 3 with ko = 2ma/Ti 2 and energy being 
in unit of koa. Inset: Solid line: 7 = 1 and other parameters 
are the same, dash line: 7 = 0.01, eVi = 30, eV2 = 0, and 
k W = 1 
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FIG. 2: Effects of ferromagnetic coupling on spin accumula- 
tion with the same parameters used in Fig. 1. Here J is in 
unit meV. 
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